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ABSTRACT 



Aims. I have analyzed a sample of seven nearby edge-on galaxies observed in the V and K'-band, in order to infer the properties of 
the dust distribution. 

Methods. A radiative transfer model, including scattering, have been used to decompose each image into a stellar disk, a bulge, and a 
dust disk. The parameters describing the distributions have been obtained through standard minimization techniques. 
Results. The dust disks fitted to the V-band images are consistent with previous work in literature: the radial scalelength of dust is 
larger than that for stars (h^/h^ ~ 1.5); the dust disk has a smaller vertical scalelength than the stellar (zi/Zs ~ 1/3); the dust disk is 
almost transparent when seen face-on (central, face-on, optical depth tq = 0.5 - 1.5). Faster radiative transfer models which neglect 
scattering can produce equivalent fits, with changes in the derived parameters within the accuracy of full fits including scattering. 
In the K'-band, no trace is found of a second, massive, dust disk which has been invoked to explain observations of dust emission 
in the submillimeter. 1 discuss the model degeneracies and the effect of complex structures on the fitted distributions. In particular, 
most bulges in the sample show a box/peanuts morphology with large residuals; two lower-inclination galaxies show a dust ring 
distribution, which could be the cause for the large fitted dust scalelengths. 
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1. Introduction 

Edge-on spiral galaxies offer a unique opportunity to study 
the three dimensional structure of disks: the vertical and ra- 
dial behavior of the stellar distribution can be analyzed in 
a rather direct manner, provided the heavily extinguished 
dust lane is avoided and/or observations are taken in a band 
not st r ongly affected by dust (s e e, e.g . [van der Kruit & Searle 
11982; de Griis&van der KruiV |199^ IPohlen etal.1 l2000bi: 
Iplorido et al. 2001); the properties of the dust distribution can 
be inferred from the extinction lane itself. In this second case, 
however, the analysis requires three-dimensional models for the 
radiative transfer of starlight through dust, which are computa- 
tionally demanding, especially when proper geometries and scat- 
tering are taken in to account. 

lb date, onlv iKvlafis & Bahcalll (1 19871) and IXilouris et all 

( 119971 119981 [19991) have fitted the surface brightness distribu- 
tion in edge-on spiral gala xies by using r ealistic radiative trans- 
fer models. In the works of lXilouris et al.L two independent disks 
were used to describe the stellar and dust distribution (both 
having a radial and vertical exponential fall-off), to which a 
spheroidal distribution was added to model the bulge. After an- 
alyzing a sample of seven galaxies, Xilouris et al. (1999) con- 
cluded that the dust disk is thinner (vertically) but larger (radi- 



ally) than the stellar; its optical depth perpendicular to the disk is 
smaller than one in optical wavebands, making the disk almost 
transparent when seen face-on. 



* Based on observations made with the Italian Telescopio Nazionale 
Galileo (TNG) operated on the island of La Palma by the Fundacin 
Galileo Galilei of the INAF (Istituto Nazionale di Astrofisica) at the 
Spanish Observatorio del Roque de los Muchachos of the Instituto de 
Astrofisica de Canarias 

** Figures |5] to [16] are only available in electronic form via 
|http://www.edpsciences.org| 



The dust disk emerging from the work of IXilouris et al.l 
poses a problem, since such a structure cannot absorb more than 
10% of the stellar radiation. Observations in the infrared, in- 
stead, show that ~3Q% of the tot al bolometric emission of a 
spiral galaxy is radiated by dust ( Popesc u & Tuffs! 12002 ) and 
thus have been absorbed from starlight. The conservation of this 
energy balance demands a larger am ount of dust in th e disk 
than what derived from optical images dBianchi et alj200b ). The 
deficit in absorption is mitigated if a sizable fraction of starlight 
is assumed to suffer extinction within localized areas (i.e. star 
forming region in molecular clouds) and not only diffuse ex- 
tinction within the disk. Even taking this into account, models 
fail to predict the emission from the colder dust component in 
the FIR/sub-mm and a supplem ental dust d istribution is needed 
(iPopescu et al.ll2000t:lMisiriotis et al.ll2001h . 

In the models fo r the edge-ons NGC 891 a nd NGC 5907, 
iPopescu et al.l (l2000l) and iMisiriotis et alj (1200 lb include a sec- 
ond dust disk, associated with the younger stellar population. 
Being thinner than the disk resulting from the analysis of 
Xilouris et al., the second disk is not detectable in optical im- 
ages, as its effects occurs in a region akeady heavily affected by 
the extinction from the first, thicker disk. Nevertheless, the sec- 
ond dust disk is required to be a factor 2-3 more massive than 
the first. An alternative explanation for the model failures in the 
FIR/sub-mm is that the adopted dust emissivity, derived from the 
diffuse e mission in the Milk y Way, is underestimated by about a 
factor 3 (iDasvra et al.ll2005h . 
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Table 1. The sample 



name 




D'' 


PA" 


RA & Dec (J2000)^ 


NGC 4013 


3 


14.5 


244.8 


11 


58 


31.3 


+43 


56 


50.8 


NGC 4217 


3 


16.9 


49.3 


12 


15 


51.1 


+47 


05 


29.8 


NGC 4302 


5 


20.2 


50.8 


12 


21 


42.3 


+ 14 


35 


53.9 


NGC 5529 


5 


41.9 


294.1 


14 


15 


34.2 


+36 


13 


37.5 


NGC 5746 


3 


26.6 


350.2 


14 


44 


56.1 


+01 


57 


18.8 


NGC 5965 


3 


47.4 


52.6 


15 


34 


02.4 


+56 


41 


07.5 


UGC 4277 


6 


76.5 


109.5 


08 


13 


57.1 


+52 


38 


53.0 



" morphology parameter (RC3). 

* distance in Mpc, assuming Hg = 73 km s"' Mpc"' and velocities 
w. r. t. the CMB dipole (from NED). 

Position angle in degrees, from the V-band fit. 
Galactic centre, from the V-band fit. 



with an offset of 100" from the galactic center. Each field 
was observed with a dithered pattern of 12 on and off posi- 
tions, each consisting of 6 short exposure of 20 sec on the 
source (and on the sky) for a total exposure of 1440 second on 
source. Due to its smaller extent on the sky, only UGC 4277 
required a single telescope pointing, with the galactic plane 
aligned with the camera diagonal. Data reduction was carried 
out with the dedicated software SNAP (Speedy Near-infrared 
data Automatic Pipeline; Mannucci et al., in preparation) which 
takes care of flat-fielding, sky subtraction, corrections for geo- 
metrical distortion and electronic effects, and final image mo- 
saicing. Calibration was obtained fr om objects in the 2MASS 
point source catalog dSkrutskie et al.ll2006) . The sky noise in the 
K'-band images is typically 21.2 mag arcsec"^ (l-o"), about a 
magnitude d eeper than in 2MASS Large Galaxy Atlas images 
(iJarrett et al.l 12003) . The seeing during observations was I'.'O 
(FWHM). 



To test the reliability of radiative transfer fitting techniques 
and ascertain if an extra dust disk component is indeed present, 
I analyze in this paper a sample of seven edge-on galaxies, ob- 
served in the V and in the K' bands. The second dust disk should 
be discernible in NIR observations, due to its higher op tical 
depth with respect to the thicker disk: iDasvra et al.l (l2005l) an- 
alyzed a K„-band image of NGC 891 and found an extinction 
lane consistent with a single disk model. In this work I will be 
able to carry the same analysis on a higher resolution and deeper 
set of NIR images. The galaxies were selected because of the 
presence of an evident dust lane, their symmetrical appearance 
and their closeness, so as to provide optimal ground for the ap- 
plication of radiative transfer fitting techniques. They all are of 
morphological type Sb to Sc. The sample is listed in Tab.[T] Two 
objects ( NGC4013 and NGC5529) have been extensively ana- 
lyzed in IXilouris et al.l (Il999t) : this will allow to check on the 
uniqueness of the results of fitting procedures. 

The paper is structured as follows: the observations are de- 
scribed in Sect.|2j the radiative transfer model and the procedure 
adopted for fitting are presented in Sect. [3j the results of fits to 
observations are shown in Sect. |4l finally, the main points are 
summarized and discussed in Sect. |5] 



2. Observations & data reduction 

Observations were carried out in March and April 2006, at the 
3.5 m TNG (Telescopio Nazionale Galileo) telescope, located at 
the Roque de Los Muchachos Observatory in La Palma, Canary 
Islands. 

V-band images were obtained in dark time using th e 
DOLORES instrument in imaging mode dMolinari et al.|[l997l) . 
The pixel scale is 0'.'275 and the field of view 9'A x 9'A, a 
size well suited to the galaxies in the sample. For each object, 
three dithered exposures were taken, for a total of 750 seconds. 
Standard data reductio n was carried out u sing the STARLINK 
package CCDPACK (.Draper et al.ll2002l) . As conditions dur- 
ing the observing nights were not photometric, calibration was 
achieved using the total V-band magnitude of the galaxies from 
RC3 dde Vaucouleurs et al. 1991). The sky noise in the V-band 
images is typically 25.8 mag arcsec"^ (l-o") and the seeing 1'.'8 
(FWHM). 

K'-band images were obtained wit h the NICS instrument in 
wide-field imaging dBaffa et al.ll2001h . Because of the smaller 
field of view (4^2 x 4'.2, with pixel scale 0'.'25), the camera 
was aligned with the galactic plane of each galaxy using the 
RC3 position angle; two overlapping fields were observed, each 



3. IVIodelling & fitting 

The analysis presented in this work consists in producing a mock 
image of a galaxy and comparing it with the observed image by 
means of minimization techniques, with the aim of deriving 
the parameters that better describe the object. 



3.1. The galactic model 

A standard description is adopted for the disk (see, e.g. 
Xilouris et al. 1999), an exponential (both along the radial co- 
ordinate r and the vertical coordinate z) with luminosity density 
(per unit solid angle) 



^disk 

p'"V,z) = :^exp 
2 Zs 



r 



(1) 



where /y'*'' is the disk /oce-on surface brightness through the 
galactic centre and hs and Zs are the radial and vertical scale- 
length. 

For the bulge I have used a de-projected Sersic profile of 
index «, with luminosity density given by 



fbulge 



bulge 



exp(-fe„ B'l") 



B" 



(2) 



An Re b/a 

where /q"'^'' is the bulge face-on surface brightness through the 
galactic centre, Rg is the effective radius, b/a is the minor/major 
axis ratio of the bulge. 



B = 



^Jr^+zy{b/a) 
R. 



/?„ = 2n - 1 /3 - 0.009876/n, a = (2n - l)/2n 

and A n ^ 2.73,3.70,4.47,5.1 2 for n = 1,2,3,4, respec- 
tively (iPrugniel & Simieiil 11997). When pr ojected on the sky 
plane, Eq. |2] corresponds to a ISersid (11968 ^ profile of index n. 
Traditio nally, bulges in spiral g alaxies have been described with 
the R^^"^ (Ide Vaucouleurs! 19591) profile typical of elliptical galax - 
ies (n - 4). Recent work, however (see, e.g. iHunt et aLll2004l) . 
suggests that bulges in late type galaxies follow preferentially 
Sersic profiles with smaller n. Here I use n - 2 and 4. 

The model includes a single dust disk with extinction coeffi- 
cient given by 



K(r, z) 



To 

2 Zd 



exp 



hd Zd 



(3) 
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where tq is the central, face-on, optical depth of the dust disk 
and /id and Zd are the radial and vertical scalelength. 

For computational reasons, the stellar disk is generally 
truncated along the radius at 4/1, and the dust disk at 4/id. 
Observations suggest that stellar disks follow the simple expo- 
nential decline of Eq.[T]up to about 4/is, although the truncation 
is not sharp and there is a large scatter in the measures (see the 
review in Pohlen et al. 2004). No indications are available for 
the dust disk. The stellar and dust disks are truncated vertically 
at 6zs and at 6zd, while the bulge extends to IQRg. 

The simulated image is created with the same pixel resolu- 
tion and extent as the observed image with which it is compared. 
The geometrical parameters describing the appearance of the 
model galaxy on the image are 0, the inclination of the galactic 
z-axis with respect to the line of sight (l.o.s; 9 - 90° for the pure 
edge-on case), the position of the projection of the galactic cen- 
ter on the sky, and the position angle (PA). PA is defined so that, 
by rotating the galaxy counterclockwise by (9Q-PA) degrees, the 
projection on the sky of the positive z-axis of the galaxy lies 
along the positive y-axis on the image. 

I have used two radiative transfer models to produce sim- 
ulated images. In the first, faster, model scattering is not taken 
into account (the extinction coefficient in Eq.[3]is taken to be en- 
tirely due to absorption). The image surface brightness is derived 
by analytic integration of the dust attenuated stellar luminosity 
pe"'" along the l.o.s. passing through the center of each image 
pixel (with the optical depth t the integral of Eq.[3]from the lo- 
cation of stellar emission to the observer). 

The second model is a Monte Carlo (MC) radiative trans- 
fer code including scattering. The basics of the method are pre- 
sented in Bianchi et al. ( 1996). The code has been rewritten and 
optimized along the lines described in Baes et al. (2003). In 
particular, the peeling-off technique tYusef-Zadeh et al.. ,1984h 
has been implemented, allowing to produce images at a spe- 
cific inclination rather than for a broad inclination band as 
in our original paper. As for the dust scattering properties, the 
albedo w and asymmetry param eter g have been taken f r om th e 
Milky Way dust grain model of IWeingartner & Draind (1200 ll) : 
for the V-band it is w = 0.67 and g = 0.54, for the K'-band 
oj - 0.45 and g = 0.14. The extinction law gives Ak' /Ay - 0. 12. 
The iHenvev & GreensteinI (fl941.) phase function for scattering 
is used. 

Finally, the model image is convolved with a gaussian of ap- 
propriate FWHM to simulate the effect of seeing. This is impor- 
tant to reduce the central spike of the bulge, especially for the 
case of the analytic calculation where, contrary to MC images, 
the surface brightness is not smoothed over the pixel area. 

3.2. The fitting procedure 

Before proceeding to the fitting, stars and image defects have 
been masked out from observations, together with all pixels lay- 
ing out of the 3-cr isophote. The sky level, derived from regions 
free of sources, has been subtracted. To speed up the fitting pro- 
cess, images are further smoothed and rebinned over 2x2 or 3x3 
pixel, typically leaving 50000 and 25000 pixels available for the 
fit, in the V- and K'-band, respectively. Tests show that this fur- 
ther rebinning does not aff'ect the fit. 

For fits with the analytic no-scattering model, I have tested 
both the Levenberg-Marquardt and the dow nhill simplex (th e 
amoeba algorithm) minimization techniques (iPress et al.lll996h . 
The first is quicker, but only when the initial guesses for the pa- 
rameters are close to the final values, otherwise the technique 
tends to be trapped into local minima of the complex surface. 



The amoeba algorithm, instead, is able, after successive restarts 
of the procedure, to converge over minima quite distant, in pa- 
rameter space, from the initial guess. When using the MC model, 
the amoeba method is mandatory, because the statistic nature of 
MC images does not allow to quickly compute the derivatives re- 
quired by the Levenberg-Marquardt method. In all cases, the fit- 
ting codes have been implemented on a parallel machine: for the 
analytic method, calculations on groups of pixels are distributed 
among different cpus and the results collected in the final im- 
age; for the MC method, each cpu produces a whole image start- 
ing from an independent random seed, and all images are finally 
summed up to produce a higher signal-to-noise result. 

To resume, the parameters obtained through the fit are 

twelve: I^''^, h„ z„ /q"'*^", ^e, b/a, hi, Zd, tq, 9, PA and the 
centre position. First, their values have been estimated from 
cu ts parallel and perpen dicular to the disk plane, as described 
in IXilouris et al.l (Il997h . Then, a no scattering fit was done, 
with successive calls of the amoeba and Levenberg-Marquardt 
methods, until a stable minimum was found. The final fit was 
achieved with the MC method, again with various calls of the 
amoeba method. 



3.3. Testing fhe procedure 

The fitting techniques were extensively tested on a large set 
of simulated images, covering a wide range of parameters. 
Simulated images were produced using the full MC model in- 
cluding dust scattering. The surface brightness in the models 
was scaled to the typical values in observations. The pixel size 
was also chosen to match the observations, and appropriate noise 
and seeing were added. Finally, the simulated images were pre- 
pared for fitting as described in Sect. 13.21 For illustrative pur- 
poses, I discuss here the results of tests on a model based on the 
NGC 891 fit of IXilouris eTall d 19991) . Simulations were made 
for the the V- and K'-band (the dust disk fitted to the V-band 
image was used in both cases, with Tq scaled according to the 
Milky Way extinction law) and for inclinations of 89.8° (the fit- 
ted galaxy inclination), 88° and 86°. 

In Fig. [1] I show the V-band surface brightness profiles of 
models (crosses) and fits (lines) for the three inclination cases. 
In each case, the profiles are obtained along three vertical cuts 
(perpendicular to the galactic plane) passing through the model 
centre and at distances from the model of and Ih^- As out- 
lined in Sect. 13.21 fits were first produced with the no-scattering 
analytical model (dashed line), then the results used as starting 
point for the MC fit (solid line). The MC fits reproduce very well 
the data and most of the input parameters can be retrieved with 
accuracies within 10%, gradually worsening for images of lower 
6 (whose value, even in the worst cases, is retrieved with an error 
of 2%). The values of and /^"'^"^ suffer the largest error (up to 
20% and 40%, respectively, in the image with = 86°). 

Fig. [T] also shows that the data can be reproduced equally 
well by the faster analytical fits which neglect scattering. The no 
scattering fit can reproduce the surface brightness of a model in- 
cluding scattering with slight changes in the parameters, which 
are still derived with the same accuracy discussed before (apart 
from Zd, /y'^'' and /q"'"*^, which can have errors as high as 20%, 
20% and 50%, respectively). The optical depth tq is underesti- 
mated (for a given optical depth, neglecting scattering results in 
an overestimate of the attenuation; Baes & Dejonghe 2001), but 
no more than 20%. In fact, th e effects of scatteri ng are minimal 
in edge-on galaxies (see, e.g. iBianchi et al.lll996ft . 
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'v.rticul o-i,c; [zj '..,:lcul o«.ct [zj 'vertical offset [z,,]' 

Fig. 1. V-band surface brightness profiles perpendicular to the galactic disk for the model and the fits described in Sect. 13.31 Three 
cuts are shown, at distances (brightest profile), 1 and 2 (dimmest profile) from the center along the galaxy's major axis. 




vertical offset [zJ vertical affset [z^] vertical offset [z^] 

Fig. 2. K'-band surface brightness profiles along the minor axis for the models with a single dust disk and the model with two dust 
disks, together with their fits. 



Models in K'-band (see Fig. |2]i have reduced extinction and 
the presence of dust can only be inferred from small asymmetries 
along the minor axis direction. Therefore, it is more difficult to 
obtain the parameters of the dust disk: the fitted /id, Zd and tq 
are within 20% of the input values. For the same reason, the 
parameters defining the stellar distributions are retrieved with 
sUghtly better accuracy than for the V-band case. 

Other systematic effects were also explored using these sim- 
ulations. Since only pixels with high S/N were chosen, a constant 
sky was subtracted from the images without attempting to fit its 
level. However, even assuming that the sky has been miscalcu- 
lated by +3cr, the parameters resulting from the fits are still very 
close to the input values, with the uncertainties quoted above. 
The same is true if the radial truncation of the stellar disk is var- 
ied within the estimates currently a vailable from observ ations. I 
confirm the correlation observed bv lPohlen et al.l (l2000al) : an un- 
derestimate in the truncation of the stellar disk causes an overes- 
timate of /is. Asymmetries and dishomogeneities in real images 
may hamper the correct determination of the inclination angle, 
especially for objects very close to edge-on: an overestimate of 
(so that the fit is more edge-on than the galaxy) would cause an 
overestimate of Zs and Zd, leading also to larger /id and smaller tq. 
It is worth noting that in the cases discussed here, where the pa- 
rameters are deliberately kept fixed on values different from the 
input ones, the residuals in the fit are generally smaller than 20% 
for most (90%) of the data points. When dealing with real im- 
ages, the deviations from the smooth model adopted here cause 
residuals of the same order or larger (Sect. |4]i. A large model 
degeneracy is thus possible. 

Finally, I have added to the simulated images of NGC 8 91 a 
second dust disk, following the recipes of lTuff^s et al.l (l2004t) . As 
a derivation of the parameters describing the second disk is not 



feasible, I have studied its effects on the single dust disk fitting 
model. The second, thin, dust disk has a central face-on, optical 
depth which is about 2.5 times that of the thick disk (tq =0.22 
in the thin disk vs tq =0.09 in the thick disk). As pointed out by 
jPopescu et alj|2000h . its effects would not be easily discerned 
from those of the thick disk in the edge-on case. However, the 
second disk would be easily seen in all cases where the first pro- 
duces a lower extinction, i.e. at lower and larger wavelength. 
Indeed, a deeper absorption trough is visible in V-band simu- 
lated images, at all inclinations. It is also visible in K'-band sim- 
ulations, as a deeper trough within 2° from the edge-on case, 
and as a stronger minor axis asymmetry at lower inclination (see 
Fig.|2]i. As Fig.|2]shows, low-residual fits can be achieved in the 
K'-band, even if the simulated image include a second dust disk: 
the parameters retrieved for the dust disk are now close to those 
of the second, more opaque, disk. 



4. Results 

For each galaxy in the sample, I show in Fig. [3] and |5] to [TOl the 
image, the fit and its residual for the V and K'-band, respectively. 
Images have been processed as described in Sect. |3.2| and rotated 
counterclockwise by 90-PA. Cuts perpendicular and parallel to 
the galactic plane, through the center and at different positions 
on the images, are shown in Fig.fflandfTTItofTSl for both bands. 
The best-fit parameters describing the stellar disk, the bulge and 
the dust disk are listed in Table |2] together with the galaxy's 
inclination 0. The galactic centre and the PA fitted to V and K' 
images generally agree within 1" and 0?5, respectively: those 
derived for V images are shown in Table [1] 
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V— band observation 




V— band model 




K' — band observation 




K' — band model 




Fig. 3. NGC 4013. From top to bottom: the V-band image, the fit to the V-band image and the residuals; the K'-band image, the 
fit to theJC-band image and the residuals. Images and models are displayed with the modified logarithmic visualization method of 
Ijarrett et al.l ( l2003i) . The scale of residuals is linear: the bar with gray shades helps to distinguish regions with values in the ranges 
[-50%,-30%], [-30%,-10%], [-10%,10%], [10%,30%] and [30%,50%] from Hghter to darker, respectively. All regions with values 
outside these ranges are coded black. The horizontal size of the box is 6'. 



6 



Simone Bianchi: The dust distribution in edge-on galaxies 




1 1 1 1 1 1 1 

■ 


1 1 1 1 1 
■ ■ ■ 


1 1 1 1 1 1 1 1 




J** 

— ^+ + 












v\ - 


- / ++ 
/ 




\ " 

-k 












































1 , , , , 


1 . . . 





vertical offset [kpc] 



radial offset [kpc] 





vertical offset [kpc] 



radial offset [kpc] 



Fig. 4. NGC 4013. Surface brightness profiles for cuts perpendicular (left panels) and parallel (right panels) to the galactic disk, for 
the V-band (top panels) and the K'-band (bottom panels) images. In the left panels, profiles are shown for distances from the fitted 
center along the galactic plane corresponding to (minor axis), +1/8 and +1/4 of the horizontal extent of the fitted area. Median 
values have been taken on a cut width of 2% that extent. The topmost data (crosses) refer to the minor axis, whose maximum has 
been assigned magnitude zero. Crosses refer to profiles on the right side of the images (with respect to Fig.O and are progressively 
displaced downward of 1.5 magnitudes when moving from the center outward. Diamonds refer to the left side profiles, and are 
further displaced by 0.5 magnitude from their specular profile on the right side. The same scheme is adopted for the surface 
brightness profiles parallel to the major axis, shown in the right panels: profiles are shown for cuts at (major axis; crosses), 1/8 and 
1/4 of the vertical extent of the fitted area, above (crosses) and below (diamonds) the galactic plane. Median values have been taken 
over a cut width of 4% that extent. To reduce the cluttering of symbols only one point every three is plotted in horizontal profiles. 
Profiles are displaced in the same manner as in the left panels. In all panels, the curves refer to analogous cuts on the fitted images. 
Gaps in both model and data correspond to masked areas (i.e. stars). The shaded boxes in the upper part of the left (right) panels 
show the position and extent of the cuts displayed in the right (left) panels. 



4.7. Global trends 

The asymmetries and dishomogeneities in the stellar and dust 
distributions result in large residuals, which also highlight the 
presence of structures different from those adopted in Sect. 13.11 
From a qualitative analysis of the fitted images and of the pro- 
files, it appears that in four objects (NGC 4013, NGC 4217, 
NGC 5529 and UGC 4277) the model can describe reasonably 
well the V images, both for what concern the stellar components 
and the extinction lanes produced by the dust disk. In two ob- 
jects (NGC 5746 and NGC 5965) the fit is marginal, with a good 
description of the general appearance of the stellar distributions, 
but with extinction features that are not reproduced by the model. 
In one object, NGC 4302, the fit is not able to describe the stellar 
distribution out of the galactic plane. In the K'-band, the general 
quality of the fits is similar to that in the V-band. 



In a fit, the parameter that most quickly converge to the final 
value is the radial scalelength /zj. In Table |2] the V-band is 
given in kpc. To allow a ready comparison, all other scalelengths 
in Table |2] have been normalized to the V-band h^. The V-band 
hs goes from about 3 kpc in the two bulge-dominated galaxies 
NGC 4013 and NGC 4217 to about 12 kpc in UGC 4277. Indeed, 
such a wide range for /is has been found when fitting the surface 
brightness distribution in large samples of less inclined galaxies 
(Ide Jondll996aHMacArthur et al.ll2003h . 

The K'-band /is is generally smaller than in the V-band, with 
a median ratio 0.75. This is in agreement with several works in 
literature, which suggest that the reduction in observed in low 
inclination galaxies when going from optical to NIR images is 
due to an intrinsic color gradie nt in the stellar population and not 
to the effects of dust extinction (Ide Jongll996bl;lMacArthur et al.l 
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Table 2. Fit results. I^''^ and are in mag arcsec , the V-band is in kpc, the K'-band and all values for Zs, Re, hd and Zd are 
in units of the V-band /zs. The Sersic index n is not fitted. The bulge-to-total light ratio B/T is computed from the stellar disk and 
bulge fit. The values for hd, Zd and 6 are derived from the V-band fit and kept fixed in the K'-band fit. 
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120031: lMollenhoffll2004l: ICunowl l2004. A notable exception is 
that oflMollenhoffet al] (l2006h . which ascribe the increase of 
radial scalelengths from B to I to the optically thick second dust 
disk of (Popescu et al. (2000). There is no evidence, instead, of 
a systematic difference between the vertical scalelength Zs in the 
V and K' -bands. In the V-band, it i s Zs/^ s ~ 0.1, a result similar 
to what obtained bv lXilouris etan(ll999h . 

More delicate is the derivation of bulge parameters. For most 
object, I have used a spheroidal bulge with Sersic index n - 2, 
a choice appropriate for the mo rpholo gical type of the galax- 
ies considered here dHunt etalJ 120041) . However, for two ob- 
jects (NGC 4013 and NGC 4217) the n = 2 bulge was too 
steep to fit both the central peak and the ellipsoidal appearance 
of more external isophotes: n - A was used in these cases. A 
spheroidal bulge appears to be a poor approximation to the real 
stellar distribution: in all but one object, positive residuals (i.e. 
model smaller than data) show 'X'-like features in the fit to 
the K' or both bands. These features reveal families of orbits 
in which disk material is t rapped by the perturbations caused 
by a bar dPatsis & Xilourisli2006), whos e edge-on view causes 
a boxy/peanut appearance dBureau et al.l 120061) . The deviations 
from the assumed bulge model are likely to affect also the deriva- 
tion of Zs (because bulge and vertical disk structure contributes 
together to the surface brightness profile out of the disk plane 
near the centre) and of the parameters describing the dust disk 
(because the effects of dust extinction are more evident in the 
central part of the disk where the bulge contribution is sizable, if 
not dominant). 

Excluding the bad, almost no-bulge, NGC 4302 fit, the n - 2 
bulges have B/T ratios in the range 0.2-0.4 in the V-band, while 
the bulge contribution is larger in the two n - A fits. In the K'- 
band fit, the B/T ratio is always smaller than in the V-band, thus 
implying that the color of the bulges is bluer than the color of the 
disks. This is puzzling, being in contrast to what generally found 
in th e bulge/disk decomposition of less inclined objects (see, e.g. 
lMollenhoff»2004i) . However, a similar variation of the B/T ratios 



can be inferred from the fit of edge-on galaxies in lXilouris et al.l 
CHI). 

The dust disks fitted to the V-band images appear to be con- 
sistent with those derived by Xilouris et al. ( 1999), although the 
spread in values is larger. Generally, the radial scalelength hd 
is larger than the stellar (Jidlh^ ~ 1.5), while the vertical scale- 
length Zd is smaller than Zs (Zs/zd ~ 3). This last condition is 
necessary to produce the well defined extinction lanes along the 
galactic disk. The disk opacity is moderate, with a V-band tq in 
the range 0.5 -1.5, making the gal axies almost transparent when 
seen face-on dXilouris et al]|l999l) . 

The extinction features in K' images are not prominent, and 
result mostly in weak asymmetries along the minor axis and in 
clumpy structures. The asymmetries due to dust cannot be easily 
distinguished by those due to deviations from the adopted mod- 
els. In most cases, a free fit would converge to a dust free model. 
Thus, the fits have been achieved by allowing tq only to vary, 
while fixing both dust disk scalelengths, hd and Zd, and the incli- 
nation 6 to the result of the V-band fit. Small values for tq, in the 
range 0.02-0.15 have been found. Such opacities only produces 
minimal extinction effects on the model (see, e.g. the solid lines 
in Fig.|2]i. 

Instead, no extinction lanes associated to a thin, massive, 
dust disk are found. If the extra dust component supposed by 
(iPopescu et al.ll20dol) exists, it does not appear to be in the form 
of a second, smooth, dust disk, which would produce more eas- 
ily discernible features (see, e.g., the dashed lines in Fig.|2]i. The 
result is in a greement with the a nalysis of a K„-band image of 
NGC 891 bv lDasvraetaIld2005l) . 

4.2. Comments on individual ga/ax/ef] 

NGC 4013: the fit is generally good for both bands, with resid- 
uals smaller than 30% in 85% of the fitted data points. The dust 
lane can be easily seen in both bands. It shows a clumpy structure 
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which lead to larger residuals along the major axis, especially in 
the V-band. A bulge with Sersic index n - A was used for this 
o bject. The V-band fit has parameters similar to those obtained 
bv lXilouris et alj (11999) . with the exception of /id, which I find 
to be smaller than hs, rather then larger This also lead to a larger 
Tq obtained here. 

NGC 4217: the same considerations as for NGC 4013 apply to 
this object. The geometrical parameters of the stell ar and dust 
disk, together with 6, are given for the B-band by lAlton et all 
(120001) : the fit here is more edge-on and is larger by 30%, 
while /zs is consistent with their fit. The difference in 9 is also 
the cause of the smaller Zs I obtain. For this object I also have 
Zd > z^'- since the bulge dominates the surface brightness apart 
from regions close to the galactic plane, an extinction lane can 
appear even if the dust disk is thicker that the stellar. The major 
axis profiles clearly show a steepening of the profile beyond 10 
kp c: this is th e signature of a Type II truncation (in the notation 
of iPohlen & Truiillo. 2006) . which cannot be reproduced by the 
sharp cut used in the model. A truncation at the s ame distance is 
also found in the Kj-band bv lFlorido et alj (1200 10 . although the 
value of /is they obtain is larger by 40% than in the present fit. 
Dust extinction can be seen in the K'-band image, although not 
in the form of a dust lane but rather as a y-axis asymmetry in 
images, because of the inclination smaller than edge-on. 
NGC 4302: this is the worst object to fit, because of the ma- 
jor axis asymmetry around the center, the presence of a sur- 
face brightness contribution from the nearby face on galaxy 
NGC 4298 (whose centre is at about 2.'3 from the galactic plane 
in the direction of the negative y axis in Fig.|6]) and because of the 
complex vertical structure. Indeed the model does not represent 
well the observed surface brightness, especially in the V-band. 
The contribution of the fitted bulge is negligible. It appears that 
there is a second, thicker stellar disk (of V-band vertical scale- 
length 1.5 kpc, a quarter of the fitted h^), with radial scalelength 
increasing with the distance from the galactic plane. Extinction 
in K' appears clumpy and does not alter significantly the sym- 
metry above the plane. 

NGC 5529: apart from a slight warp in the external part of the 
disk, a fit can be achieved, with residuals smaller than 30% in 
70% of the data points. The fit of the disk is consisten t with 
the analysis of a V-band image in IXilouris et al] (Il999l) . Only 
tenuous, clumpy extinction can be seen in the K' image. 
NGC 5746: the fit quality is the same as for NGC 5529, in both 
bands. The fit to the extinction lane in the V-band, though, is 
poorer The parameters for the stellar and du st disks are not con - 
sistent with those reported for the B-band by I Alton et alj (l2000h . 
with /is larger by 40% in the present case. Extinction is located 
in a clumpy ring structure of radius a!l.'3, which is evident in the 
K' image. From the ring in K', an inclination 6 - 86° + 1° can 
be estimated, which is consistent with the fitted value. 
NGC 5965: the galaxy is warped, with the outer disk and a 
ring-like dust lane on a different plane with respect to the inner 
peanut-shaped bulge. The stellar ring is evident in the K' im- 
age, while no extinction features can be detected. As with NGC 
5746, the fit to the extinction lane in the V-band is poor, despite 
the global fit quality is not much different from the other ob- 
jects. An estimate for /is is available for theJCj-band which is 
about 20% smaller than the present value dBizvaev & Mitronoval 
12002 ). The fitted 6 is consistent with the aspect of the K' ring, 
which suggests = 83.5° + 1.5°. 

UGC 4277: the fit is good, with 90% of the data points fitted by 
the model with residuals smaller than 30%. While the dust lane 
is evident for this nearly edge-on object in the V-band, in the K'- 
band only a weak asymmetry can be seen for the bulge above 



an d below the galactic plane. Th e K-^band value for hs derived 
by ' Bizvaev & Mitronoval (l2002h is consistent with the result of 
the present fit. 



5. Summary & discussion 

In the present work a sample of seven nearby edge-on galax- 
ies, observed in the V and K'-band, has been analyzed using a 
radiative transfer model, in order to constrain the parameters de- 
scribing the stellar and dust distributions. The dust and stellar 
disks have been fitted with two smooth, independent, exponen- 
tial disks, and the bulge with a spheroidal distribution following 
Sersic profiles with n = 2 and 4. Hereafter I summarize and dis- 
cuss the main results of the work: 

(i) When applying the fitting technique to mock images of the 
same extent and S/N as the observations, it became evident 
that the parameters cannot be constrained better than a few 
tens percent, because of the model degeneracies. In particu- 
lar, it is possible to obtain good fits even using faster radiative 
transfer model which neglect scattering. For V-band mod- 
els, this leads only to a small underestimate of the dust disk 
opacity, and to variation of the other fitted parameters within 
the accuracy of full fits including scattering. Furthermore, 
the differences between scattering and non scattering mod- 
els are smaller than the typical residuals obtained in fits of 
real images. 

(ii) In the V-band, the parameters obtained for the dust disk 
are generally consistent with those derived bv lXilouris et al.l 
(.1999.) on a different sample of objects: the dust disk has a 
larger radial scalelength than stars (hd/hs ~ 1.5); it has a 
smaller vertical scalelength (zd/Zs ~ 1 /3); it is almost trans- 
parent when seen face-on (tq = 0.5 - 1.5). However, discrep- 
ancies exists when fits to single objects are compared to other 
works in literature, possibly a result of different models, fit- 
ting techniques and image coverage. It is worth noting that 
the largest values for have been obtained for two objects, 
NGC 5746 and NGC 5965, which clearly show the presence 
of a ring structure in the K'-band image. If dust is distributed 
in a ring (indeed, extinction along the ring is seen in the K' 
image of NGC 5746), the displacement of the maximum ex- 
tinction out of the center of the galaxy could be the reason 
for the larger /id, when the fit is performed using the simple 
exponential disk of Eqn. |3] 

A similar effect could be caused by the more external clouds, 
if the clumpy struct ure of the dust disk domina tes over its 
smooth component. Misiriotis & Bianchi (2002) have ana- 
lyzed the effects of clumping on radiative transfer fits and 
concluded that a clumpy distribution leads to a smaller h^. 
The different result may be due to the fact, however, then 
their model was unable to produce high resolution images 
with discernible dust clouds, 
(iii) The properties of the stellar disk are similar to those derived 
from bulge/disk decomposition of less inclined galaxies. In 
particular, the stellar disk has a larger in the V-band than in 
the K'-band, suggesting that the color gradient along disks is 
due to the intrinsic properties of the stellar populations rather 
than to the effect of extinction. The bulge component proved 
to be more complex than the assumed spheroidal model, with 
an 'X'-like residual betraying a box/peanut morphology in 
all but one galaxy. I caveat here that the departure of the 
bulge shape from model assumptions can cause uncertain- 
ties in the derivation of the dust disk properties, since the ex- 



Simone Bianchi: The dust distribution in edge-on galaxies 



9 



tinction lane is deeper in regions where the bulge contributes 
most. 

(iv) Images in the K'-band show little dust extinction, mostly 
consisting in minor axis asymmetries and in clumpy struc- 
tures. I have found no trace of the second, thinner, more 
ma ssive and smooth du st dis k which has been intr oduced 
by Popes cu et"an ([2000') and 'Misiriotis et alj (1200 ll) to ex- 
plain the sub-mm dust emission in NGC 891 and NGC 5907. 
Since the problem of the discrepancy in energy balance be- 
tween the radiation absorbed by dust from stars and that 
emitted in the infrared persists, w e are left with the second 
hypothesis in Popescu et al.l ('2000): dust could be preferen- 
tially distributed in optically thick quiescent or star-forming 
clouds. Furthermore, dust in dense clouds could have a larger 
emissivity than what usually assumed for Milky Way dust, 
thus reheving part of t he energy balance discrepancy (see, 
e.g. iDasvra et al.ll2005l for an alternative solution to the de- 
ficiency of sub-mm emission in models of edge-on galax- 
ies). Realistic models for dust emission in a complex disk 
are needed, in order to ascertain the relative contribution of 
difiiise and clumpy dust to the infrared spectral energy dis- 
tribution. 

Clearly, the uncertainties in the derivation of the structure 
of the dust disk could be reduced by using a dust tracer more 
direct than its extinction on starlight. In particular, observa- 
tions of dust emission in the sub-mm appear to be more easy 
to interpret in terms of dust density, because of the relatively 
high resolution and of the small dependence of the Rayleigh- 
Jeans spectrum on the dust temperature gradient. In the next 
few years, the advent of high-sensitivit y, second generatio n sub- 
mm bolomet er-arrays U ke SCU BA-2 ([Holland etal.1 1200 6) and 
LABOCA ( Krevsa etani2003l) . will allow to perform observa- 
tions of spirals over a much wider extent along the disk than 
previously possible. Equally important are the infrared observa- 
tions of nearby galaxies carried on with the Spitzer satellite: a 
large number of nearby spirals has bee n already observed for the 
SINGS survey dKennicutt et al.ll2003h . and analysis is ongoing 
on a sample of local edge-on objects (Roelof de Jong and Benne 
Holwerda, private communication). Very promising for a high 
resolution study of the dust distribution appears to be the corre- 
latio n between the PAHs emissio n, the ISM density (iRegan et al.l 
I2OO61) and the sub-mm emission dHaas et al.ll2002l) . 
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5. NGC 4217. Image order is the same as in Fig. [3] The horizontal size of the box is 8' 
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6. NGC 4302. Image order is the same as in Fig. [3] The horizontal size of the box is 8'. 
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7. NGC 5529. Image order is the same as in Fig. [3] The horizontal size of the box is 7'. 



Simone Bianchi: The dust distribution in edge-on galaxies, Online Material p 5 




Fig, 8. NGC 5746. Image order is the same as in Fig. [3] The horizontal size of the box is 8'. 
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. 9. NGC 5965. Image order is the same as in Fig. [3] The horizontal size of the box is 6'. 
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10. UGC 4277. Image order is the same as in Fig. [3] The horizontal size of the box is 4'. 
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11. Same as Fig.g] but for NGC 4217. 
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12. Same as Fig.g] but for NGC 4302. 
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Fig. 14. Same as Fig.g] but for NGC 5746. 
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15. Same as Fig.g] but for NGC 5965. 
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16. Same as Fig.g] but for UGC 4277. 



